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The calculation of the cosmological relic density of the dark matter candidate within supersymmetric models is an 
interesting possibility to obtain additional constraints on the supersymmetric parameter space with respect to col- 
lider, electroweak precision, and low-energy data. Considering future cosmological precision measurements, radiative 
corrections can play an important role in the analysis. We present full QCD and SUSY-QCD corrections to neutralino 
pair annihilation into quark-antiquark pairs and analyze their impact on the neutralino annihilation cross section and 
the relic density. By comparing to the relic density obtained from a pure leading order calculation, we show that the 
corrections strongly influence the extraction of SUSY mass parameters from cosmological data. 



1. INTRODUCTION 



Beginning with the rotational curves of galaxies, the astronomical and astrophysical evidence gathered in recent 
(-H ' years supports the hypothesis of dark matter in our Universe. Moreover, the models of structure formation favor 
04 the so-called cold dark matter which consists of weakly interacting particles with non-relativistic velocities. Based 
^ • mainly on the cosmic microwave background measurements, one obtains the following limits on the relic density of 
O , cold dark matter cross-section 

-2 



0.095 < 17cdm/i < 0.136, (1) 

where the 2a interval was obtained from the three- year data of the Wilkinson Microwave Anisotropy Probe (WMAP), 
combined with the results of Sloan Digital Sky Survey (SDSS), Supernovae Legacy Survey (SNLS) experiment and 
from baryonic acoustic oscillations data [if. 

One of the main motivations for physics beyond the Standard Model (SM) is the inability of SM to provide candidates 
for cold dark matter. In the Minimal Supersymmetric Standard Model (MSSM) with R-parity conservation, the 
lightest supersymmetric particle (LSP) is always stable. If it is also a colour singlet, electrically neutral and massive, 
] it is a good candidate for cold dark matter. In many scenarios, in particular those where supersymmetry breaking 
^ is gravity mediated, the lightest neutralino is the LSP and a suitable dark matter candidate. We can then calculate 
r '"j its number density n by solving the Boltzmann equation 

where the first term on the right-hand side corresponds to a dilution due to the expansion of the Universe, and the 
second one to a decrease due to annihilations and co-annihilations. Here, H denotes the Hubble expansion parameter, 
7icq the density of the relic particle in thermal equilibrium and v is the relative velocity of the annihilating pair. 
The Boltzmann equation can be integrated numerically and the energy density of the dark matter particle is approx- 
imately inversely proportional to the thermally averaged annihilation cross-section (CTann^^)i 

^CBMh^ oc r . (3) 

The latter includes all annihilation and co-annihilation processes of the dark matter particle into Standard Model 
particles. These by now standard calculations are incorporated in program packages such as DarkSUSY 2] and 
MicrOMEGAs which allow different scenarios and different models to be studied. 

Using the experimental limits ([T]), we can put constraints on the MSSM parameter space. If one is confined to the 
MSSM with gravity mediated supersymmetry breaking (mSUGRA) , the constraints allow for a rather narrow allowed 
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mo (GeV) 


mi/2 (GeV) Ao 


(GeV) tan/3 sgn(/i) 




bb 


tt 


1 


1500 


130 


1500 


10 
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0.116 


86% 
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1500 


1500 





50 


+ 


0.112 


83% 
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5300 


625 


1500 


10 


+ 


0.110 




72% 


4 


3000 


600 





50 
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0.110 


15% 


64% 



Table I: The minimal supergravity (mSUGRA) scenarios that lead to a neutralino relic density agreeing with the current limits 
of Eq. (Ill, satisfy the current SUSY-particle mass limits, and present important contributions of neutralino annihilation into 
bottom or top quark-antiquark pairs. We indicate the mSUGRA parameters at the high scale, the relic density obtained with 
micrOMEGAs, and the contributions of the different quark final states to the annihilation cross-section {(Tannt'). 

band in the space of the high-scale parameters mg, mi/2, ^Oj tan/3, sgn//. Moreover, the constraints are expected to 
become even more strict in the future when the Planck mission will provide additional data. 

To match and fully exploit the current and expected experimental precision, one has to improve the theoretical 
prediction of the relic density by including the radiative corrections. We calculated the QCD and SUSY-QCD 
corrections to neutralino annihilation into quarks. In the work presented here, we investigate the effects of radiative 
corrections on the constraints of the MSSM parameter space. 

2. NEUTRALINO ANNIHILATION INTO QUARKS 

Neutralino annihilation into fermions is one of the most important annihilation channels. It is dominant for light 
neutralinos where all other final states are not kinematically allowed. The leading order cross-section was discussed 
in much detail in The cross-section in the non-relativistic limit is conventionally written as (see 0]) 

va^a + bv^ + 0{v'^) , (4) 

where v is the relative velocity of the neutralino pair. In the case of annihilation into fermions, the coefficient a is 
proportional to the mass of the fcrmion 4]. As a consequence, the dominant final states are the third generation 
fermions, i.e. in our case the top and bottom quarks. Note that the contribution from the coefficient b is suppressed 
with respect to the one coming from a. 

The annihilation into fermions has to be enhanced by a resonance in order to yield a small enough relic density 
satisfying the experimental constraints ([T]). It proceeds mainly through a Higgs boson resonance (in particular the 
CP-odd Higgs resonance). This implies that Yukawa coupling of fermions to the Higgs boson becomes an essential 
parameter in this calculation along with tan/3. The bottom Yukawa coupling receives large SUSY corrections for 
large tan/? which is another reason to incorporate radiative corrections in the predictions. The approximation 
used in I6( is already included in the MicrOMEGAs package alongside with the QCD corrections for Higgs decays into 
fermions (see Q). 

3. RESULTS 

We calculate the complete QCD and SUSY-QCD corrections not only to the Higgs boson exchange but also to 
the s-channel Z-boson exchange and the t-channel and u-channel squark exchange. We also include parts of the 
correction for the Higgs exchange which are not treated in [il . We find that the additional correction when including 
the full one-loop radiative result are 10-15% relative to the corrections already included in MicrOMEGAs. 
We present numerical results for four typical mSUGRA parameter points shown in Tab.[ll The high-scale parameters 
were evolved down to the electroweak scale using SPheno and the neutralino relic density has been calculated with 
micrOMEGAs, where we have included our calculation of the annihilation cross section as discussed above. For each 
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m^=1500, m^^=!30. A^=-150Q. lar^=10, n>0 m^^lSOO GeV. A ^-1500 GeV. iah^^W,\i>0 




p^^(GeV) m,„(GeV) m„ (GeV) 

Figure 1: The effects of radiative corrections on cross-section as a function of centre-of-mass momentum (left) and the effects 
on the resulting dark matter's relic density as a function of mSUGRA parameters mi/2 (middle), mo (right). The results are 
shown for parameter points 1 and 2 with a dominant annihilation into bottom quarks. 

point, the relic density fulfills the limits ((T|) and additionally, the branching ratio of the decay b ^ sj lies within the 
experimental limits 

BR(6-^S7) = (3.55 ± 0.26) • 10~^ (5) 

obtained from combined measurements from BaBar, Belle, and CLEO This limit leads to constraints on param- 
eters connected to sfermions and gauginos. 

The selected points with a dominant bottom quark final state are either near the bulk region (point 1) with low 
fermion mass parameter mi/2 for low tan/3 (the scalar mass is rather large in order to avoid the constraint by 6 — » 57) 
or in the Higgs-funnel region for large tan /3 where our point 2 is positioned exactly on the CP-odd Higgs resonance. 
Both points 3 and 4 for which the top quark final state is important, lie in the focus point region where the scalar 
mass parameter mo is large. 

In Figs. [T] and Owe show the numerical results for the relevant annihilation cross section and the resulting neutralino 
relic density for the points presented in Tab. HI We show the leading order result (green dash-dotted), the approx- 
imation already implemented in micrOMEGAs (red dashed), and the result including our full one-loop SUSY-QCD 
corrections (blue solid) . In the graphs showing the cross section (left row) we also indicate the Boltzmann distribution 
function involved in the calculation of the thermal average (couloured area). 

As already stated above, the one-loop contributions that are not implemented in micrOMEGAs lead to an increase of 
about 10 - 15 % in the annihilation cross section both for bottom and top quark final states. In consequence, the 
neutralino relic density decreases by about the same amount for the different points, as can be seen in the center 
and right panels showing the neutralino relic density as a function of the high-scale mass parameters mi/2 and mo, 
respectively. For fixed values of tan/3, Aq and sgn/z, the favoured band in the mo-m]^y2 plane is shifted to lower 
gaugino masses and higher scalar masses (focus point region). 
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Figure 2: The effects of radiative corrections on cross-section as a function of centre-of-mass momentum (left) and the effects 
on the resulting dark matter's relic density as a function of mSUGRA parameters mi/2 (middle), mo (right). The results are 
shown for parameter points 3 and 4 with a dominant annihilation into top quarks. 
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